Cosmic vector for dark energy: constraints from SN, CMB and BAO 
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It has been recently shown that the presence of a vector field over cosmological scales could 
explain the observed accelerated expansion of the universe without introducing neither new scales 
nor unnatural initial conditions in the early universe, thus avoiding the coincidence problem. Here, 
we present a detailed analysis of the constraints imposed by SNIa, CMB and BAO data on the vector 
dark energy model with general spatial curvature. We find that contrary to standard cosmology, 
CMB data excludes a fiat universe for this model and, in fact, predicts a closed geometry for the 
spatial sections. We see that CMB and SNIa Gold data are perfectly compatible at the 1-sigma 
level, however SNIa Union dataset exhibits a 3-sigma tension with CMB. The same level of tension 
is also found between SNIa and BAO measurements. 
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I. INTRODUCTION 

More than ten years have passed since the first indi- 
cations of the accelerated expansion of the universe |lj 
and still today it remains as one of the most intriguing 
problems in cosmology. Moreover, the accelerated ex- 
pansion has been confirmed during the last decade by 
many different probes, mainly through measurements of 
Type la supernovae, CMB temperature power spectrum 
and the Baryon Acoustic Oscillations (BAO) @, |3j 
Since Einstein's gravity predicts that a universe contain- 
ing (baryonic and dark) matter and radiation should be 
decelerating rather than accelerating, these observations 
could be signalling cither the breakdown of General Rela- 
tivity at cosmological scales or the existence of some sort 
of non-ordinary energy with negative pressure known as 
dark energy. 

The simplest model describing dark energy is the exis- 
tence of a cosmological constant. Although it fits obser- 
vations with very good precision, it suffers from theoret- 
ical problems since the inferred scale for the cosmologi- 
cal constant turns out to be around 10~'^eV, which is a 
very tiny value compared to the gravitational scale set by 
Newton's constant G ~ Mp"^ with Mp - lO^^GeV. This 
poses a problem of naturalness because of the existence 
of two scales in the theory which differ in many orders of 
magnitude. Moreover, a related problem arises because 
the amount of energy density stored in the form of dark 
energy is comparable to that stored in the form of mat- 
ter at the present epoch, in spite of having evolved very 
differently in the past. Thus, in order to get these prob- 
lems around several models have been proposed which are 
mostly based on either cosmological scalar fields 0, jc 1 or 
infrared modifications of the gravitational action [7, Q 
(see Q for a review on dark energy models). However, 
none of these models succeeded in solving the previously 
addressed problems of the cosmological constant because 
they either introduce new dimensional scales in the action 
or unnatural initial conditions to get the right accelera- 



tion at the right time. 

In a previous paper [lo| we showed that a model based 
on the dynamics of a vector field on cosmological scales 
can give rise to a period of accelerated expansion without 
introducing neither new dimensional scales nor unnatu- 
ral initial conditions, thus avoiding fine-tuning or coin- 
cidence problems. Although the existence of periods of 
accelerated expansion for vector field models was already 
known the model proposed in [l^l consists just of 

the simplest kinetic terms for the vector field containing 
two fields and two derivatives and without any potential 
for the field. In fact, in this paper we shall show that 
such a model is nothing but a gauge field with a gauge 
fixing term coupled to the Ricci tensor. Moreover, very 
recently [l^ it has been also shown that the own elec- 
tromagnetic field could be a natural candidate for dark 
energy. 

In order to constrain dark energy models from obser- 
vations, we typically use distance indicators so that we 
can confront distance measurements to the corresponding 
model predictions fT3|. To do that we can resort to two 
different types of objects, namely: standard candles and 
standard rulers. Standard candles are objects of known 
intrinsic luminosity, so that the corresponding comoving 
distance can be determined. That way, it is possible to 
reconstruct the Hubble expansion rate by searching this 
sort of objects at different redshifts. The most important 
class of such indicators are Type la supernovae. On the 
other hand, standard rulers are objects whose comoving 
size is known, so that we can measure the angular dis- 
tance and, therefore, compare to that predicted by the 
dark energy model. A well-known example is the sound 
horizon size at the last scattering surface. This standard 
ruler can be measured directly from the CMB temper- 
ature power spectrum and, also, from Baryon Acoustic 
Oscillations (BAO) through the matter power spectrum 
at low redshift. 

The paper is organized as follows: in section I we in- 
troduce the vector model for dark energy and derive all 
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the necessary equations for the rest of the paper. In sec- 
tion II we show how the different distance indicators wiU 
be used in order to obtain constraints for the model. Fi- 
nally, section III is devoted to the results obtained from 
the analysis. 



II. VECTOR DARK ENERGY 

The action proposed in [l^] to describe dark energy is 
the following: 



S = / (Tx^/^ - 



R 



- ^V^A^V^A'' + ]^R^,A^^A''^ (1) 

However, this action can be written in a more suggestive 
form by integrating by parts and taking into account that 
R^,^A^'A'' = (V^A'')2 - V^AyVA^'. Then, the action 
for the vector field can be expressed as: 



S 



1 -f^M"-^ 



R 



i {V^^A^^f + R^.A'^A'^ 



(2) 



In this form, it becomes apparent that the theory is that 
of a gauge vector field with a gauge-fixing term (in the 
Feynman gauge) and coupled to the Ricci tensor. This 
coupling provides an effective mass term driven by grav- 
ity. 

Now, we take variations with respect to the metric and 
the vector field to obtain the Einstein's and vector field 
equations respectively: 



R^v 



1 
2 

R^..A'' = 



(3) 

(4) 



where T^i, is the conserved energy-momentum tensor for 
matter and radiation (and all other possible components) 
and is the energy-momentum tensor for the vector 
field. We consider the case in which the vector field is 
homogeneous and only has time-component, i.e., A^ — 
(A)(i), 0, 0, 0). However, unlike in the previous paper 
[Kll, here we will consider the effects of the curvature so 
the metric is given by: 



dt^ - a{tf 



dr^ 



1 — fcr^ 



(5) 



In this metric, equation ([4]) for the homogeneous time- 
component of the vector field reads: 



Ao + iHAo - 3{2H^ + H] Aq ^ 



(6) 



where H = d/a is the Hubble parameter. For our pur- 
poses on this work it will be useful to express the latter 



equation in terms of the redshift 



1 as follows: 



d^An 



dz^ (1 + z)H{z) dz 



H{z) 



[l + zf 



;|l(l + =)'.4ol^0 



(7) 



This equation can be easily solved when the universe is 
dominated by radiation or matter, assuming that the con- 
tribution of the vector field is negligible. In those epochs 
the Hubble parameter is given by H = p/t ~ p{\ -l-z)^/^ 
with p = 1/2 for radiation and p = 2/3 for matter and, 
according to ([7]), the vector field evolves as: 



A^{z) = A+{l + zr++A^{l + zr- 



(8) 



with Aq constants of integration and a± = (1 ± l)/2 in 
the radiation era, and a± = (3 ± -v/sS) /4 in the matter 
era. Notice that, since the vector field is constant during 
the radiation era, we can set the initial conditions at any 
time during that epoch without modifying the evolution 
of the universe. 

On one hand, the o) component of Einstein's equa- 
tions can be written as: 

-2 = n,rr{l+zf + nr{l+zf + nk{l + zf + pA{z) (9) 
^0 

where is the present value of the Hubble parameter, 
which is usually expressed as Hq = 100 km s~^Mpc~^, 
VLm and ilr are the density parameters corresponding to 
matter and radiation respectively, Vlk — and pA is 

the energy density of the vector field, whose expression 
is: 



PA = 



IP 

Hi 



\aI -ii 



dAg 

dz 



(10) 

Moreover, we measure the vector field in units of the 
reduced Planck mass Mp = I/vSttG. We can obtain 
the evolution of this energy density in the radiation and 
matter dominated eras by using the growing solution in 
([8]) for the vector field to give: 



PA = PAo{l + zY 



(11) 



where k = 4 in the radiation era and k = (9 — v33) /2 ~ 
1.63 in the matter era. Notice that the energy density 
scales as radiation in the early universe so that pa/ Pb, 
is constant in that epoch (see in Fig. [1]). Therefore, 
we obtain once again that the time at which we set the 
initial conditions in the early universe for the vector field 
lacks importance, since its fraction of energy density is 
constant during that epoch. 

Finally, when the vector field becomes dominant the 
Universe suffers a Type HI singularity Ts*! in which its 
evolution ends at a finite time t^nd but with a finite size 
a^nd- Moreover, as we approach the final singularity we 
have pDE — *■ oo and poE — cxd, with wde — oo, 
whereas the vector field takes a finite value. 
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FIG. 1: Left panel: evolution of energy densities. Dashed (red) for radiation, dotted (green) for matter and solid (blue) for 
vector dark energy. We also show for comparison the cosmological constant energy density in dashed-dotted line. We see the 
rapid growth of dark energy contribution at late times approaching the final singularity. Right panel: cosmological evolution 
of the vector field. In these two plots we see the unimportance of the time at which we set the initial conditions due to the fact 
that both the vector field and the fraction of dark energy density are constant in the early universe. 



On the other hand, the (' i) component of Einstein's 
equations is: 

H"^ \ 1 + z dH^ , -.A , .0 

= nr{i + zY + nk{i + zf-PA{z) 



i/2 



3- 



l + zdH 



(12) 

where we have used that pr = \pr and p,„ = and the 
pressure of the vector field is given by: 



Pa{z) 



^0 



5 4. \ dH 



42 
^0 



, , ^ dAo 3 , , n / dAo 



.(13) 



In order to perform the analysis in next sections it will 
be convenient to write equations ^ and in terms of 
{rirh^ , ^Imh^ , ^kh^} as follows: 



H 



n^h'^{i+zf+nrh'^{i+z)^+nkh'^{i+z)^+PA{z) 

(14) 



l + zdH 
H dz 



nrh^{i+z)^+nkh'^{i+zf-PA{z) 

(15) 

where H = i7/(100km s-^Mpc-^), i.e., H{z = 0) = /i. 
Note that neither pA nor pA depend on the normalization 
of the Hubble parameter. Moreover, flrh^ contains the 
contribution of photons as long as relativistic neutrinos. 



rirh^ ^ n^h^il + 0.2271iVe//) 



(16) 



with Neff = 3.04 the effective number of neutrino species 
and n-yh^ = 2.469 x 10"^ for Tcmb = 2.725. 



The model is completely determined once wc fix the 
set of parameters {$7™, il^, A^ad} so that the model has 
three free parameters. To confront the model to SN and 
BAO dataset we only need to integrate the system of 
equations up to redshift ~ 2 whereas CMB dataset re- 
quires to obtain the solution up to the last scattering 
surface so that the method to solve the equations will be 
different for each case. The present value of the Hub- 
ble expansion rate is no longer a free parameter in this 
model because it can be obtained in terms of the previ- 
ous parameters after integrating the equaions. In fact, we 
could take {rirm^kT h} as independent parameters and, 
therefore, Arad would already be determined, although 
this approach is more difficult to implement numerically. 
Notice that this model contains exactly the same number 
of parameters as ACDM. 



III. LIKELIHOOD CALCULATIONS 

In this section we shall explain the procedure followed 
to confront the vector dark energy model to the different 
distance indicators. 



SN 



The apparent magnitude of a supernova placed at a 
given redshift z is related to the expansion history of the 
universe through the distance modulus: 



p : 



M — 51ogDi — 51og/i -|- po 



(17) 
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where m and M are the apparent and absolute magni- 
tudes respectively, /ig = 42.38 and Dl = HqcIl with cIl 
the luminosity distance = (1 + z)r(z), being r(z) the 
comoving distance, given for the metric ^ by: 



iz) = 



Hp 
H{z') 



dz' 



(18) 



with Sk[x] — sinx, x, sinhx for Q^. < 0,flk = 0,17^ > 
respectively. 

Then, to confront the model to each supernovae data 
set we construct the corresponding estimator: 



X^SN 



1=1 



(19) 



which must be marginalized over h in order to obtain the 
constraints on the parameters Qm and Qk- 

In order to calculate Xsn' ^^"^ ^^'^^ ^^at the 

SNIa dataset corresponds to redshifts below 2 so that we 
can neglect the contribution from radiation in Einstein's 
equations. With this in mind, we solve numerically the 
system of equations ([7]) and (fT2|) for H/Hq and Aq. As 
this system is of second order with respect to Aq and 
first order with respect to H/Hq we need to set the intial 
values of Ag, dA^/dz and H/Hq. However, these three 
initial values are related by means of Friedmann equation 
© so that we can obtain the initial value for H/Hq in 
terms of the initial values of Aq and its derivative. On 
the other hand, as we know the analytic solution of the 
vector field in the matter dominated era as that given 
in ([8]) , we can relate the initial value of the derivative 
of the vector field to the initial value of the vector field 
(neglecting the decaying mode). Therefore, we only need 
to give the initial value for Aq in order to set the initial 
conditions and we are left with A™^, flm and ilfe as free 
parameters in terms of which we obtain the correspond- 
ing XsN estimator. Therefore, we shall use: 



ifJ.{zi; fir, 

i=l 



^k,Ar 



(20) 



instead of (jl9p and marginalize over Aini. 

In this work we have used two sets of supernovae: the 
Gold set j2] and the more recent Union set Q . 



B. BAO 

BAO measurements provide the following distance ra- 
tios [H: 



Vbao ^ ( -Z^^ 



0.1980 ± 0.0058 
^_ , , 0.1094 ±0.0033 

Dv'(0.35) 

where rs{z) is the sound horizon size given by: 
1 /■T+^ da 

rs(z) 



, (21) 



(22) 



and 



Dv{z)= rHz) 



H 



1/3 



(23) 



is the dilation scale. Finally, Zd is the drag epoch at 
which baryons were released from photons and which can 
be calculated by using the fitting formula 20]: 



Zd 

with 
bi 



2\0.251 



1 + 0.659(r2„/i2) 



0.828 



[i + bi{nbh 



2\b2 



(24) 



0.313(a„/i2 



-0.419 



[1 + 0.607(f}„/i2) 



2^0.6741 



52 = 0.238(f7„/i2)0-223. 
Then, we define the BAO array: 



Dv(0.2) 
Dv(0.35) 



- 1.980 
-0.1094 



so that: 



Xbao — X^J^oC^^qXbao- 



In this expression, the inverse covariance matrix is 



p-i 

^BAO 



35059 -24031 
-24031 108300 



(25) 
(26) 



(27) 



(28) 



(29) 



The procedure we follow in this case is analogous to that 
used for the SNIa analysis, although, as Xbao depends 
on the amount of baryons Qbi we also need to marginalize 
over this parameter. 



CMB 



Following [17[, we use the distance priors method to 
confront dark energy models to CMB data [ll|. This 
method uses two distance ratios measured by means of 
the CMB temperature power spectrum: 

• The "acoustic scale", which measures the ratio of 
the angular diameter distance to the decoupling 
epoch and the comoving sound horizon size at de- 
coupling epoch. This first distance ratio can be 
expressed as: 



7rr(z*) 
rs{z^) 



(30) 



Moreover, we use the fitting formula of proposed 
in fii: 



1048 [1 + 0.00124(176/1^) -"-738] 
X [l + gi(r!„/i2)92] 



with 



91 



92 



0.0783(f7f,/i2 



-0.238 



1 + 39.5(J7f,/l2)0.763 

0.560 

' l + 21.1(17fcft2)1.81 



(31) 

(32) 
(33) 
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• The second distance ratio measures the ratio of the 
angular diameter distance and the Hubble radio at 
the decoupling time. It is usually called the "shift 
parameter" and can be expressed as 

R - ^n„,H^r{z,) (34) 

The values reported in fl7t| for these distance priors 
are: 

/ Ia{z*) \ ( 302.10 ±0.86 \ 

McMB= R{z*) = 1.710 ±0.019 (35) 
\ j \ 1090.04 ±0.93 / 

with the following inverse of the covariance matrix 

/ 1.800 27.968 -1.103 \ 
^cliB = 27.968 5667.577 -92.263 . (36) 
\ -1.103 -92.263 2.923 / 

Then, we define the CMB array as 

( Ia~ 302.10 \ 
XcAfB= i?- 1.710 , (37) 
\z^- 1090.04 ) 

so that: 

x\;mb —^^MB^CMB^cm-B- (38) 

The procedure we follow in this case is somewhat dif- 
ferent to that used in the previous sections. The main 
difference comes from the fact that CMB distance pri- 
ors are evaluated at a time when radiation is important 
so that we cannot neglect its contribution in Einstein's 
equations anymore. To simplify numerical calculations 
we use equations (O, (fT4|) and (fT5|) . Notice that, un- 
like the SN and BAO approach, from these equations 
we obtain the Hubble expansion rate normalized to 100 
km s~'^Mpc~^ so that H{z = 0) = /i. Thus, for given 
{Vl„ih? ,^kh^} we use ([T4|) to relate the initial condition 
for the Hubble expansion rate to the intial condition of 
the vector field (that we shall name Arad) and, then, 
solve numerically ([7]) and (|15p. Since the initial condi- 
tions are set in the radiation-dominated era when, ac- 
cording to (IHI), the vector field is constant, the initial 
condition for the derivative of the vector field is set to 
zero. Moreover, the constancy of the vector field dur- 
ing that epoch eliminates the dependency on the time 
at which we place the initial conditions, i.e., Arad does 
not depend on Zini. That way, we obtain the expansion 
rate H{z) that will allow us to compute the distance in- 
dicators described above in terms of {Q,mh^ , ^kh^, Arad} 
(notice that such indicators do not depend on the nor- 
malization of the Hubble expansion rate). Hence, we can 
compute the corresponding Xcmb which will depend on 
{rim^i^, ^fcfe^, ^bh^, Arad} and, following the prescription 
given in [l7[, we marginalize over Qbh^ and Arad (which 



is equivalent to marginalize over h) and use the result- 
ing marginalized likelihood to obtain the corresponding 
contours. 

Since CMB distance priors were derived in [l7j assum- 
ing that dark energy was not important at decoupling 
time (z* ~ 1090) and given that the vector field model 
does not produce a significant amount of dark energy at 
high redshifts, these priors are, in principle, applicable in 
this case. 



IV. RESULTS 

In this section we present the results obtained after 
confronting the model given by the action ^ with the 
tests explained above. We have also performed the anal- 
ysis for a ACDM model for comparison. 

Using the Gold data set we obtain a best fit for flm = 
0.385 and flA = 0.611 with xLn = 172.92, which is the 
same value found in where we imposed flat spatial 
sections. This is understandable because, from the above 
values of and we obtain flk = 0.0043 so that the 
best fit is very close to the flat case. However the la 
contour allows both open and close universes and, unlike 
(lol | , a wide range of values for f2m and Qa is within the 
1(7 region, as we can see in Fig. For a ACDM model 
with non-vanishing curvature we obtain the best fit for 
f]„ = 0.46 and = 0.98 with xliin = 175.04 so we 
still obtain a better fit to the Gold data set than ACDM. 
On the other hand, the best fit obtained for the vector 
dark energy model from the Union data set corresponds 
to n,n 0.260 and ^a = 0.503 with xLn = 311.96. 
From Fig. ^ we see that this data set favors an open 
universe for this model, being the flat case at more than 
2cr. For ACDM the best fit happens for = 0.41 and 
f^A = 0.93, being Xmm = 310.23, which is lower than 
that obtained for the vector field. This effect is probably 
due to the SNLS points contained in the Union data set 
which, as it is shown in [13], favor ACDM over the vector 
field model at more than 2a in the fiat case. However, 
when the fiatness assumption is dropped, ACDM fits the 
Union data set better than the vector field model only at 
less than la. 

Concerning BAO dataset, it favors an open universe 
with a small amount of matter for the vector field model, 
as we see in Fig. ©. Moreover, the compatibility of these 
data with SNIa data is only at the 3a level. However, it is 
worth mentioning that these distance indicators are ob- 
tained after analyzing the actual observational data with 
ACDM as fiducial model so that its applicability to test 
dark energy models is justified as long as such models do 
not differ much from a cosmological constant. Neverthe- 
less, this is not the case for the vector dark energy model 
whose equation of state varies very rapidly and, indeed, 
has a future singularity so that the obtained 3a tension 
could be due to the dependency of BAO data on the fidu- 
cial model. In any case, this is the less confident dataset 
to constraint the vector model and, in general, any dark 
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FIG. 2: In this plots we show the 68%, 95% and 99% C.L. regions for BAO (orange), CMB (green) and SNIa (blue). We show 
the contours obtained for both the Union data set (left) and the Gold data set (right). The blue line corresponds to a flat 



universe. 




FIG. 3: In this plot we show the likelihood obtained from the 
CMB dataset for the vector field model. We can see that a 
flat universe is clearly ruled out and a closed geometry for the 
spatial sections is strongly favored. 



energy model, since it may give shifted parameters due 
to a biased determination of the sound horizon scale due 
to the presence of additional relativistic degrees of free- 
dom, early dark energy or a non-standard recombination 
scheme [2l|. 

Finally, CMB data is totally incompatible with flat 
spatial sections and, in fact, it predicts a closed universe 
with a wide range of flm allowed. These results show 



that, contrary to common belief, CMB data do not nec- 
essarily favors a flat universe. In Fig. [3l the correspond- 
ing likelihood for the CMB dataset is plotted and we can 
see how the flat case is ruled out for the vector model. 

In Fig. [2] we see that CMB contours are compatible 
with BAO at 2a level for small values of flm and flk 
close to zero. Concerning SNIa contours, CMB is in con- 
flict with the Union data set contours at more than 3cr 
whereas it is compatible at la level with the Gold data 
set. 



V. CONCLUSIONS 

In this work we have performed a detailed analysis of 
the constraints imposed by SNIa, CMB and BAO data 
on the vector dark energy model proposed in p^ . We 
have considered cosmologies with arbitrary spatial cur- 
vature and obtained confidence regions in the (flrm^A) 
plane. We have found that for the SNIa Gold dataset, 
the vector model fit is better than that of ACDM, but for 
the Union dataset the situation is reversed. We find that 
contrary to standard cosmology, CMB data excludes a 
flat universe for this model and, in fact, predicts a closed 
spatial geometry. On the other hand, CMB and SNIa 
Gold data are perfectly compatible at the 1-sigma level, 
however SNIa Union dataset exhibits a 3-sigma tension 
with CMB. The same level of tension is also found be- 
tween SNIa and BAO measurements, although this may 
be due to the dependency of BAO measurements on the 
fiducial model. 
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